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ABSTRACT 

Context. ( Ori A is a hot star claimed to host a weak magnetic field, but no clear magnetic detection was obtained so far. In addition, 
it was recently shown to be a binary system composed of a 09.51 supergiant and a BHV star. 

Aims. We aim at verifying the presence of a magnetic field in ^ Ori A, identifying to which of the two binary components it belongs 
(or whether both stars are magnetic), and characterizing the field. 

Methods. Very high signal-to-noise spectropolarimetric data were obtained with Narval at the Bernard Lyot Telescope (TBL) in 
Erance. Archival HEROS, FEROS and UVES spectroscopic data were also used. The data were first disentangled to separate the two 
components. We then analyzed them with the Least-Squares Deconvolution (LSD) technique to extract the magnetic information. 
Results. We confirm that ( Ori A is magnetic. We find that the supergiant component f Ori Aa is the magnetic component: Zeeman 
signatures are observed and rotational modulation of the longitudinal magnetic field is clearly detected with a period of 6.829 d. This 
is the only magnetic O supergiant known as of today. With an oblique dipole field model of the Stokes V profiles, we show that the 
polar field strength is ~ 140 G. Because the magnetic field is weak and the stellar wind is strong, ( Ori Aa does not host a centrifugally 
supported magnetosphere. It may host a dynamical magnetosphere. Its companion ( Ori Ab does not show any magnetic signature, 
with an upper limit on the undetected field of ~ 300 G. 
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1. Introduction 


Magnetic fields play a significant role in the evolution of hot 
massive stars. However, the basic properties of the magnetic 
fields of massive stars are poorly known. About 7% of the mas¬ 
sive stars are found to be magnetic at a level that is detectable 
with current instrumentation (Wade et al. 2014| l. In particular, 
only 11 magnetic O stars are known. Detecting magnetic field 
in O stars is particularly challenging because they only have 
few, often broad, lines from which to measure the field. There is 
therefore a deficit in the knowledge of the basic magnetic prop¬ 
erties of O stars. 

We here study the O star ^ Ori A. A magnetic field seems to 
have been detected in this star by Bouret et al. (2008 1 . Their de¬ 
tailed spectroscopic study of the stellar parameters led to the de¬ 
termination of an effective temperature of - 29500 + 1000 
K and log g = 3.25 + 0.10 with solar abundances. This makes 
^ Ori A the only magnetic O supergiant. Moreover, [Bouret et al.| 
( |2008| l found a magnetic field of 61 + 10 G, which makes it 
the weakest ever reported field in a hot massive star (typically 
ten times weaker than those detected in other magnetic massive 
stars). They found a rotational period of ~1 days from the tempo¬ 
ral variability of spectral lines and the modulation of the Zeeman 
signatures. To derive the magnetic properties, they used six lines 


* Based on observations obtained at the Telescope Bernard Lyot 
(USR5026) operated by the Observatoire Midi-Pyrenees, Universite de 
Toulouse (Paul Sabatier), Centre National de la Recherche Scientifique 
of Prance. 


that are not or only weakly affected by the wind. The rotation pe¬ 
riod they obtained is compatible with their measured v sin /=100 
km s“'. 


In addition, the measurement of the magnetic field provided 
by jBouret et al. ( 2008| l allows characterizing the magnetosphere 
of ^ Ori A and locating it in the magnetic confinement-rotation 
diagram ( Petit et al. 12013) 1: ( Ori A is the only known magnetic 
massive star with a confinement parameter below 1, that is, with¬ 
out a magnetosphere. 


For all these reasons, the study of the magnetic field of 
f Ori A is of the highest importance. Each massive star that is 
detected to be magnetic moves us closer to understanding the 
stellar magnetism of hot stars. Studying this unique magnetic 
massive supergiant is also of particular relevance for our under¬ 
standing of the evolution of the magnetic field in hot stars. 


^ Ori A has a known BOIII companion, ^Ori B. In addition, 
Hummel et al. ( 2013|l found that ^ Ori A consists of two com¬ 


panion stars located at 40 mas of each other, orbiting with a pe¬ 
riod of 2687.3+7.0 days. To determine a dynamical mass of the 
components, Hummel et al.| (2013 1 analyzed archival spectra to 
measure the radial velocity variations. The conclusions reached 
are presented below. The primary ( Ori Aa is a 09.51 supergiant 
star, whose radius is estimated to 20.0+3.2 Rq and whose mass 
is estimated to 33+10 Mq. The secondary ^OriAb is a BlI'V 
with an estimated radius of 7.3+1.0 Rq and an estimated mass of 
14+3 Rq. Moreover, fOri A is situated at a distance of 387 pc. 
Initial estimates of the elements of the apparent orbit were ob¬ 
tained by [Hummel et al.|(|2013|l using the Thiele-Innes method. 
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The estimation provided a value of the periastron epoch of JD 
2452734.2+9.0 with a longitude of 24.2+1.2°. The eccentricity 
is estimated to be 0.338+0.004. 

Bouret et al. ( 2008) 1 considered ^OriA as a single star of 
40 Mq with a radius equal to 25 Rq, seen from Earth at an in¬ 
clination angle of 40°. Taking into account that the star is a bi¬ 
nary could strongly modify the magnetic field value derived for 
only one of the binary components. In their analysis, the mag¬ 
netic signature was normalized by the full intensity of the lines 
from both components, and if only one of the two stars is mag¬ 
netic, the field was thus underestimated. Moreover, the position 
in the magnetic confinement-rotation diagram will be modified 
as a result of the new magnetic strength value, but also as a con¬ 
sequence of the new stellar parameters. 

Based on new spectropolarimetric observations of fOriA 
and archival spectra presented in Sect.[2| we here seek to confirm 
that ^ Ori A is a magnetic star (Sect. 1^ We determine with sev¬ 
eral techniques, including by disentangling the composite spec¬ 
trum (Sect. 1^ whether the magnetic field is hosted by the pri¬ 
mary or the secondary star of (Ori A. We then determine the 
field strength of the magnetic component (Sect. and quan¬ 
tify the non-detection of a field in the companion (Sect. [^. In 
addition, we investigate the rotational modulation of the mag¬ 
netic field, its configuration (Sect.|^, and the possible presence 
of a magnetosphere (Sect.J^. Finally, we discuss our results and 
draw conclusions in Sect. IW 


2. Observations 


2.1. Narval spectropolarimetric observations 

Spectropolarimetric data of 4" Ori A were collected with Narval 
in the frame of the project Magnetism in Massive Stars (MimeS) 
(see e.g. Neiner et al.||2011 1. This is the same instrument with 
which the magnetic field of ^ Ori A was discovered by |Bouret| 
et al. (20081. Narval is a spectropolarimeter installed on the two- 
meter Bernard Lyot Telescope (TBL) at the summit of the Pic du 
Midi in the French Pyrenees. This fibre-fed spectropolarimeter 
(designed and optimized to detect stellar magnetic fields through 
the polarization they generate) provides complete coverage of 
the optical spectrum from 3700 to 10500 A on 40 echelle orders 
with a spectral resolution of ~65000. Considering the size of the 
fiber, the light from ^ Ori B was not recorded in the spectra, but 
light from both components of ^ Ori A was collected. 

^ Ori A was first observed in October 2007 during 7 nights 
(PI: J.-C. Bouret) and these data were used in [Bouret et al.| 
( 2008| l. Then, this star was observed again in October 2008 dur¬ 
ing 5 nights (PI: J.-C. Bouret) and between October 2011 and 
February 2012 during 24 nights by the MiMeS collaboration 
(PI: C. Neiner). This provides a total number of 36 nights of 
observations. The observations were taken in circular polarimet- 
ric mode, that is, measuring Stokes V. Each measurement was 
divided into four subexposures with a different polarimeter con¬ 
figuration. 

Since {Ori A is very bright (V=1.77), only a very short ex¬ 
posure time could be used to avoid saturation. To increase the 
total signal-to-noise ratio (S/N), we thus obtained a number of 
successive measurements each night, which were co-added. The 
exposure time of each subexposure of each measurement varies 
between 20 and 120 s, and the total integration time for a night 
varies between 1280 and 7680 s (see Table[^. 

Data were reduced at the telescope using the Libre-Esprit re¬ 
duction package ( jPonati et al.|1997] l. We then normalized each 
of the 40 echelle orders of each of the 756 spectra with the con- 


Table 1. Journal of Narval observations. The columns list the 
date and the heliocentric Julian date (HJD) for the middle of 
observation, the number of sequences, and the exposure time per 
individual subexposure, the signal-to-noise ratio in the I profiles, 
and the orbital phase. 


# Date 


1 

17oct07 

2 

18oct07 

3 

19oct07 

4 

20oct07 

5 

21oct07 

6 

23oct07 

7 

24oct07 

8 

22oct08 

9 

23oct08 

10 

24oct08 

11 

25oct08 

12 

26oct08 

13 

04octll 

14 

05octll 

15 

lOoctll 

16 

lloctll 

17 

30octll 

18 

07novll 

19 

llnovll 

20 

12novll 

21 

24novll 

22 

25novll 

23 

26novl1 

24 

29novll 

25 

30novll 

26 

14decll 

27 

08janl2 

28 

13janl2 

29 

14janl2 

30 

15janl2 

31 

16janl2 

32 

25janl2 

33 

26janl2 

34 

08febl2 

35 

09febl2 

36 

10febl2 


mid-HJD 

2454391.559 

2454392.719 

2454393.570 

2454394.491 

2454395.518 

2454397.496 

2454398.526 

2454762.644 

2454763.645 
2454764.654 
2454765.639 
2454766.635 
2455839.688 
2455840.670 
2455845.608 
2455846.632 
2455865.712 
2455873.557 
2455877.626 
2455878.565 
2455890.673 
2455891.660 
2455892.502 
2455895.667 
2455896.600 
2455910.477 
2455935.555 
2455940.536 
2455941.539 
2455942.475 
2455943.367 
2455952.529 
2455953.431 
2455966.472 
2455967.402 
2455968.343 


Texp is) 

48 X 4 X 20 
8 X 4 X 40 
44 X 4 X 40 
48 X 4 X 40 
48 X 4 X 40 
48 X 4 X 40 
48 X 4 X 40 
40 X 4 X 50 
38 X 4 X 50 

36 X 4 X 50 

37 X 4 X 50 

38 X 4 X 50 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 

5 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 

6 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
12 X 4 X 90 
8 X 4 X 90 
12 X 4 X 90 

11 x4x 120 

12 X 4 X 90 


S/N l^orb 


4750 

0.617 

2220 

0.617 

6940 

0.617 

6860 

0.618 

7070 

0.618 

7180 

0.619 

7270 

0.619 

6660 

0.755 

5530 

0.755 

6790 

0.756 

6140 

0.756 

6420 

0.756 

5810 

0.156 

5790 

0.156 

2040 

0.158 

3450 

0.158 

5610 

0.165 

2700 

0.168 

4860 

0.170 

4830 

0.170 

4180 

0.175 

4900 

0.175 

4490 

0.175 

5400 

0.176 

2030 

0.177 

1360 

0.182 

5630 

0.191 

5060 

0.193 

5350 

0.193 

4680 

0.194 

4520 

0.194 

5120 

0.198 

3200 

0.198 

3900 

0.203 

4340 

0.203 

2198 

0.203 


tinuum task of IRAlQ Finally, we co-added all the spectra ob¬ 
tained within each night to improve the S/N, which varies be¬ 
tween 1360 and 7270 in the intensity spectra (see Table [2. We 
therefore obtained 36 nightly averaged measurements. 


2.2. Archival spectroscopic observations 

In addition to the spectropolarimetric data, we used archival 
spectroscopic data of (Ori A taken with various echelle spec¬ 
trographs. 

In 1995, 1997 and 1999, spectra were obtained with the 
HEROS instrument, installed at the ESO Dutch 0.9 m telescope 
at the La Silla Observatory. The spectral resolution of HEROS is 
20000, with a spectral domain from about 350 to 870 nm. In ad¬ 
dition, in 2006, 2007 and 2009, data were taken with the EEROS 
spectrograph installed at the ESO 2.2 m at the La Silla observa¬ 
tory. The spectral resolution of EEROS is about 48000 and the 
spectral domain ranges from about 370 to 900 nm. Finally, in 

’ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 


2 



























A. Blazere et al.: The magnetic field of f Ori A 


21oct07 


23oct07 


4e-05 
2e-05 
V 0 
-2e-05 
-4e-05 
4e-05 
2e-05 
N 0 
-2e-05 
-4e-05 





22oct08 



-400 -200 0 200 400 

24oct08 


5e-05 ^ 
V 0 


-5e-05 ^ 
5e-05 


N 0 










-400 -200 ,,0,.200 400 -400 -200 ,.0,.200 400 -400 -200 ,.0 , 200 400 -400 -200 ,.0 , 200 400 


V (km/s) 


V (km/s) 


V (km/s) 


V (km/s) 


Fig. 1. LSD Stokes I (bottom), Stokes V (top), and null N (middle) profiles, normalized to Ic, for 8 selected nights. The red line is a 
smoothed profile. 


Table 2. Journal of archival spectroscopic observations of 
^OriA obtained with HEROS, EEROS and UVES, showing 
the date and heliocentric Julian date, instrument used, exposure 
time, signal-to-noise ratio, and orbital phase. 


Date 

JD 

Instrument 

Texp 

S/N 

0orb 

1995 

2449776.024 

HEROS 

57x1200 

1200 

0.90 

1997 

2450454.379 

HEROS 

16x1200 

1000 

0.15 

1999 

2451147.333 

HEROS 

64x1200 

1200 

0.41 

2006 

2453738.159 

EEROS 

60 

100 

0.37 

2007 

2454501.018 

EEROS 

2x20 

250 

0.66 

2009 

2454953.970 

EEROS 

5x10 

200 

0.84 

2010 

2455435.373 

UVES 

36x2 

2000 

0.01 


2010, spectra were taken with the UVES spectrograph ( |Dekker| 
et al. 12000 1 installed at the VET at the Paranal Observatory. Its 
spectral domain ranges from about 300 to 1100 nm with a spec¬ 
tral resolution of 80000 and 110000 in the blue and red domains 
respectively. 

We co-added spectra collected for each year to improve the 
final S/N. We therefore have seven spectra for seven different 
years, with a S/N of between about 100 and 2000 (see Table|^. 


3. Checking for the presence of a magnetic fieid 


The magnetic field of ^ Ori A claimed by |Bouret et al. ( |2008| l has 
not been confirmed by independent observations so far and one 
of the goals of this study is to confirm or disprove its existence 
using additional observations. 

To test whether (Ori A is magnetic, we applied the Least- 
Squares Deconvolution (LSD) technique (|Donati et al.||1997)l. 


We first created a line mask for ^ Ori A. We started from a list 
of lines extracted from VALD ([Piskunov et al. 1995[ Kupka & 


Ryabchikova 1999| l for an O star with Teff^JOGOO K and log 


g=3.25, with their Lande factors and theoretical line depths. We 
then cleaned this line list by removing the hydrogen lines, the 
lines that are blended with hydrogen lines, and those that are not 
visible in the spectra. We also added some lines visible in the 
spectra that were not in the original O-star mask. Altogether, we 
obtained a mask of 210 lines. We then adjusted the depth of these 
210 lines in the mask to fit the observed line depths. 

Using the final line mask, we extracted LSD Stokes 1 and V 
profiles for each night. We also extracted null (N) polarization 
profiles to check for spurious signatures. The LSD I, Stokes V 
and the null N profiles are shown in Fig.[2for 8 of 36 nights. A 
plot of all profiles is available online (Fig. [TT). Zeeman signa¬ 
tures are clearly seen for these 8 nights and some others as well, 
but are not systematically observed for all nights. We calculated 
the false alarm probability (FAP) by comparing the signal inside 
the lines with no signal ( |Donati et al. 119971 1. If FAP < 0.001%, 
the magnetic detection is definite; if it is 0.001% < FAP < 0.1% 
the detection is marginal, otherwise there is no detection. Table|^ 
indicates whether a definite detection (DD), marginal detection 
(MD) or no detection (ND) was obtained for each of the night. 
DD or MD were obtained in 15 out of 36 nights. The existence of 
Zeeman signatures confirms that ^ Ori A hosts a magnetic field, 
as previously reported by Bouret et al. (2008|l. 


The previous study of the magnetic field of ( Ori A ( [Bouret 
[et al.|[2008] l only used a few lines that were not affected by 
the wind. However, we need to use as many lines as possible 
to improve the S/N. We therefore checked whethe our results 
were changed by using lines that might be affected by the wind. 
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Fig. 2. LSD Stokes I (bottom), Stokes V (top), and null N (mid¬ 
dle) profiles, normalized to Ic, for the night of 25 October 2008 
for the mask of |Bouret et al.| ( |^08| ) (left) and for our own mask 
(right). The red line is a smoothed profile. 


We compared the LSD results obtained with the mask used by 
Bouret et al. (20081 and our own mask (see Fig. |^. The sig¬ 


natures in Stokes V are similar with both masks and the mea¬ 
surements of the longitudinal magnetic field are consistent (e.g., 
44.5+19.6 G with the mask of [Bouret et al.| ( |2008l l and 35.8+7.2 
G with our mask for measurement #11, see Fig. |^. However, 
the S/N is better with our mask (the S/N of Stokes V is 57624) 
than with the mask of jBouret et al.| ( [T008| l (the S/N of Stokes V 
is 27296). Therefore, we used all available lines for this study. 

However, the line mask used in this first analysis includes 
lines from both components of ^ Ori A. We thus ignored which 
component of the binary is magnetic or whether both compo¬ 
nents are magnetic. To provide an answer to this question, we 
must separate the composite spectra. 


4. Separating the two components 

4.1. Identifying the lines of each component 


We hrst created synthetic spectra of each component. The 
goal was to identify which lines come from the primary com¬ 
ponent, the secondary component, or both. To this aim, we 
used TLUSTY ( jHubeny & Lanz 1995| l. This program cal¬ 
culates plane-parallel, horizontally homogeneous stellar atmo¬ 
sphere models in radiative and hydrostatic equilibrium. One of 
the most important features of the program is that it allows 
for a fully consistent, non-LTE metal line blanketing. However, 
TLUSTY does not take winds into account, which can be impor¬ 
tant in massive stars, especially in supergiants. 

For the primary component ^ Ori Aa, we computed a model 
with an effective temperature 7’eff=29500 K and logg=3.25, cor¬ 
responding to the spectral type of the primary as given by 


Hummel et al. (20131. For the secondary, we computed a model 


with an effective temperature 7’eff=29000 K et logg=4.0, again 
following Hummel et al. P013| l. We used solar abundances, for 
both stars. The emergen spectrum from a given model atmo¬ 
sphere was calculated with SYNSPECj^ This program is com¬ 
plemented by the program ROTINS, which calculates the rota- 


^ Synspec is a general spectrum synthesis pro¬ 
gram developed by Ivan Hubeny and Thierry Lanz: 
http://nova.astro.umd.edu/Synspec49/synspec.html 


tional and instrumental convolutions for the net spectrum pro¬ 
duced by SYNSPEC. 

Comparing these two synthetic spectra to the observed spec¬ 
tra of Ori A, we identified which lines belong only to ^ Ori Aa, 
only to ^ Ori Ah, and which are a blend of the lines of both com¬ 
ponents. If one observed line only existed in one synthetic spec¬ 
trum, we considered that this line is only emitted from one com¬ 
ponent of the binary. If it existed in both synthetic spectra, we 
considered this line to be a blend of both components. We then 
created line lists containing lines from the three categories (only 
Aa, only Ah, or both). 

In addition, we gathered archival spectra of ^ Ori A taken 
with the spectrographs EEROS, HEROS and UVES (see 
Sect. |2.2| l. While these spectra do not include polarimetric in¬ 
formation, they cover the orbital period much better than our 
Narval data (see Tableand Fig.[^. In particular, some spectra 
were obtained close to the maximum or minimum of the radial 
velocity (RV) curve. 

We compared the spectrum taken close to the maximum and 
minimum RV, because the line shift is maximum between these 
two spectra. We arbitrarily decided to use the spectra taken close 
to the maximum as reference. Depending on the shift, we deter¬ 
mined the origin of the lines. If the lines of the spectrum taken 
at minimum RV are shifted to the blue (respectively red) side 
compared to the spectrum at maximum RV, the line comes from 
the primary Aa (respectively secondary Ab) component. When 
lines from the two components are blended, the core of the lines 
are shifted to the red side and the wings to the blue side. 

The identification of lines made this way resulted in similar 
line lists as those obtained by comparing the observed spectra 
with synthetic ones. 

We then ran LSD again on the observed Narval spectra, 
once with the mask containing the 157 lines only belonging to 
(Ori Aa and once with the mask only containing the 67 lines 
from /"OriAb. We observe magnetic signatures in the LSD V 
prohles of ^ Ori Aa that are similar to those obtained in the orig¬ 
inal LSD analysis presented in Sect. In the contrast, we do 
not observe magnetic signatures in the LSD Stokes V prohles of 
Ori Ab. We conclude that ^ Ori Aa is magnetic and ( Ori Ab is 
not. 

However, the LSD prohles of ^ Ori Aa obtained this way are 
very noisy, because of the low number of lines in the mask, and 
they cannot be used to precisely estimate the longitudinal mag¬ 
netic held strength. To go further, it is necessary to disentangle 
the spectra, so that more lines can be used. 


4.2. Spectral disentangling of Narval data 

We hrst attempted to use the Fourier-based formulation of the 
spectral disentangling (hereafter, spd) method (Hadrava 1995| l 
as implemented in the FDBinary code ( Ilijic et al.||2004 1 to si¬ 
multaneously determine the orbital elements and the individual 
spectra of the two components Aa and Ab of the Ori A binary 
system. The Fourier-based spd method is superior to the original 
formulation presented by Simon & Sturm| ( [l994| l that is applied 
in the wavelength domain in that it is less time-consuming. In 
particular, this increases the technique’s efficiency when it is ap¬ 
plied to long time-series of high-resolution spectroscopic data. 

One of the pre-conditions for the spd method to work ef- 
hciently is a homogeneous phase coverage of the orbital cy¬ 
cle with the data. In particular, covering the regions of maxi¬ 
mum/minimum radial velocity (RV) separation of the two stars 
is essential, because these phases provide key information about 
the RV semi-amplitudes of both stellar components. 
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Fig. 3. Orbital phase distribution of the spectra of ^ Ori A. The 
diamonds indicate the radial velocity measured by 
( 2013| l. The lines correspond to the best fit of the 
of each component. The crosses correspond to our Narval ob¬ 
servations and the triangles to the archival spectroscopic data. 
Phase zero corresponds to the time of periastron passage (To = 
2452734.2 HJD) as reported by |Hummel et al.| ( |2()13| l. 


Hummel et al. 
radial velocity 


Figurej^illustrates the phase distribution of our Narval spec¬ 
tra according to the orbital period of 2687.3 days reported by 
Hummel et al.| ( |20T3] l. Obviously, the spectra provide very poor 
phase coverage; no measurements exist at phases ~0.0 and 0.4, 
corresponding to a maximum RV separation of the components 
(see Fig. 5 in |Hummel et al.|2013 i. This prevents determining of 
accurate orbital elements from our Narval spectra. 

Our attempt to use the orbital solution obtained by Hummel| 
et al. ( 2013| l to separate the spectra of the individual compo¬ 
nents also failed; although all regions in which we disentangled 
the spectra indicate the presence of lines from the secondary in 
the composite spectra, the separated spectra themselves are un¬ 
reliable. 


4.3. Disentangling using the archival spectroscopic data 


Since the spd method failed in disentangling the Narval spectra 
because of the poor phase coverage, we again used the spectro¬ 
scopic archival data obtained with FEROS, HEROS, and UVES. 
The orbital coverage of these spectra is much better than the one 
of the Narval data. We have seven spectra taken at different or¬ 
bital phases, including phases of maximum RV separation of the 
components (see Table |^. We first normalized the spectra with 
IRAE. We used the orbital parameters given by [Hummel et al.| 
( 2013| l for the disentangling. 

The coverage of these spectra enables the disentangling us¬ 
ing EDBinary. As an illustration of the results, a small part of the 
disentangled spectra is shown in Eig.|^ The results confirms the 
origin of the lines that were determined in Sect 4.1 and also the 
spectral types of the components given byjHummel et al.|(|2013|l. 


5. Measuring the longitudinal magnetic field of 
Ori Aa 

Eollowing these results, we assume that ^OriAb is not mag¬ 
netic and that the Stokes V signal only comes from ^ Ori Aa. 
Therefore, we ran the LSD technique on the Narval spectra with 
a mask containing all lines emitted from ^ Ori Aa, even those that 


are blended with the ones of ^ Ori Ab, to obtain the LSD Stokes 

V profile of ^ Ori Aa. The contribution of ^ Ori Ab to this Stokes 

V signal will be null, as the magnetic signal is only provided by 
( Ori Aa. 

However, we were unable to disentangle the Narval data (see 
Sect. |4.2| l, therefore the LSD Stokes I spectra of ^Ori Aa could 
not be computed in the same way as the LSD Stokes V spectra. 
As a consequence, we attempted to compute the LSD Stokes I 
profiles in several ways. 


5.1. Using the Narval data and correcting for the companion 


Eor the I profiles, we first proceeded in the following way: we 
computed the LSD Stokes I profiles with different masks that 
only contained the lines of ^ Ori Aa, only the lines of ^ Ori Ab, 
and the only blended lines. We subtracted the LSD Stokes I pro¬ 
files obtained for the lines of ( Ori Ab alone from the LSD Stokes 
I profiles obtained for blended lines to remove the contribution 
from the Ab component. We then averaged the LSD Stokes I 
profiles obtained this way and the one obtained for the lines of 
f Ori Aa alone. In this way the same list of lines (those of Aa 
alone and the blended ones) are used in the final LSD Stokes I 
profiles as in the LSD Stokes V profiles calculated above. 

This allowed us to use more lines than in Sect. |4.1| (i.e., to 
include the blended lines) and to improve the resulting S/N. We 
obtained magnetic signatures similar to those derived in Sects. 
and 4.1 (see Eig. [^l. However, the S/N remained low, and some 
contribution from the Ab component is probably still present 
in the LSD Stokes I profile. Longitudinal field values extracted 
from these LSD profiles may thus be unreliable. 


5.2. Using synthetic intensity profiles 

To improve the LSD I profiles, we attempted to use the synthetic 
TLUSTY/SYNSPEC spectra calculated in Sect. [4.1 [ for ^ Ori Aa. 
We ran the LSD tool on the synthetic spectra to produce the syn¬ 
thetic LSD Stokes I profiles of ^ Ori Aa with the same line mask 
as the one used for the LSD Stokes V profiles above. 

We then computed the longitudinal magnetic field values 
from the observed LSD Stokes V profiles and the synthetic LSD 
Stokes I profiles. We calculated the longitudinal magnetic field 
B; for all observations with the center-of-gravity method ([Rees 
|& Semelj 1979) 1, 

f vV(v)dv 

Bi = -2.14 X 10^"- ^ - G. 

^ogmCjU - I{y))dv 

We obtained longitudinal magnetic field values between -144 
and H-112 G with error bars between 20 and 100 G. 


5.3. Using disentangled spectroscopic data 

Although we were unable to disentangle the Narval data, we 
obtained disentangled spectra from the purely spectroscopic 
archival data. To derive the longitudinal magnetic field values 
more accurately, we therefore used the disentangled spectra ob¬ 
tained from the purely spectroscopic data. 

We ran the LSD technique on the disentangled archival spec¬ 
tra obtained for Ori Aa using the same line list as we used for 
Stokes V. Thus, we obtained the observed mean intensity pro¬ 
file for Ori Aa alone. We then computed the longitudinal mag¬ 
netic field values from the observed LSD Stokes V profiles from 
Narval and the observed LSD Stokes I profiles from the disen¬ 
tangled spectroscopic spectra. 
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wavelength (angstrom) 


Fig. 4. Small part of the spectrum of ^ Ori A showing the com¬ 
posite observed spectrum (black), the spectrum of ^'OriAa 
(blue) and that of ( Ori Ah (red). 


The shape of the magnetic signatures in LSD Stokes V pro¬ 
files (Fig. 1^ is similar to the shapes obtained for the combined 
spectra (Fig. [TJ and the various methods presented above. The 
LSD Stokes I spectra now better represent the observed ^ Ori Aa 
spectrum, however. We therefore adopted these LSD profiles in 
the remainder of this work. Fifteen of the 36 measurements are 
DD or MD. 

As above, we calculated the longitudinal magnetic field B/ 


for all observations with the center-of-gravity method (Rees & 
Semel \919). Results are given in Table The longitudinal 


field B; varies between about -30 and -1-50 G, with typical error 
bars below 10 G. N values are systematically compatible with 
0 within 3crN, where crN is the error on N, while B/ values are 
above 3crB/ in seven instances, where crB/ is the error on B/. 


6. No magnetic field in ^Ori Ab 

6.1. Longitudinal magnetic field values for ( OriAb 

To confirm the non-detection of a magnetic field in ^ Ori Ab, we 
ran the LSD technique with a mask that only contained lines 
emitted from ^OriAb, that is 67 lines. This ensures that the 
LSD Stokes V profiles are not polluted by the magnetic field 
of ^Ori Aa. Signatures in the LSD Stokes V profiles are not de¬ 
tected in any of the profiles (all ND), as shown in Table|^and in 
Fig. I^for selected nights when a signal is detected in ( Ori Aa. 

Using these LSD profiles and the center-of-gravity method, 
we calculated the longitudinal field value, the null polarization, 
and their error bars for ^ Ori Ab. We find that both B/ and N 
are compatible with 0 within 3cr for all nights (see Table |^. 
However, the error bars on the longitudinal field values of 
^ Ori Ab are much higher (typically 70 G) than those for ( Ori Aa 
(typically 10 G), because far fewer lines could be used to extract 
the signal for ( Ori Ab. 

6.2. Upper limit on the non-detected field in ^ OriAb 

The signature of a weak magnetic field might have remained hid¬ 
den in the noise of the spectra of the ^ Ori Ab. To evaluate its 
maximum strength, we first fitted the LSD I profiles computed 
above for ^ OriAb with a double Gaussian profile. This fit does 
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Fig. 5. Examples of LSD Stokes I profiles (bottom) computed 
from the disentangled spectroscopic data, Stokes V (top) and 
null N (middle) profiles, normalized to Ic, from the Narval data 
for the primary component ^ Ori Aa for a few nights of observa¬ 
tions. The red line is a smoothed profile. 


not use physical stellar parameters, but it reproduces the I pro¬ 
files as well as possible. We then calculated 1000 oblique dipole 
models of each of the LSD Stokes V profiles for various values of 
the polar magnetic field strength Bpoi. Each of these models uses 
a random inclination angle i, obliquity angle f5, and rotational 
phase, as well as a white Gaussian noise with a null average 
and a variance corresponding to the S/N of each observed pro¬ 
file. Using the fitted LSD I profiles, we calculated local Stokes 

V profiles assuming the weak-field case, and we integrated over 
the visible hemisphere of the star. We obtained synthetic Stokes 

V profiles, which we normalized to the intensity of the contin¬ 
uum. These synthetic profiles have the same mean Lande factor 
and wavelength as the observations. 

We then computed the probability of detecting a dipolar 
oblique magnetic field in this set of models by applying the 
Neyman-Pearson likelihood ratio test (see e.g. Helstrom]|1995[ 
Kay|199^ |Levy|2008| l. This allowed us to decide between two 


hypotheses: the profile only contains noise, or it contains a noisy 
Stokes V signal. This rule selects the hypothesis that maximizes 
the detection probability while ensuring that the EAP is not 
higher than 10“^ for a marginal magnetic detection. We then cal- 
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Table 3. Longitudinal magnetic field of the magnetic primary 
star ^ Ori Aa. The columns list the heliocentric Julian dates 
(HJD) for the middle of observation, the longitudinal magnetic 
field and its error in gauss, the detection status: definite detection 
(DD), marginal detection (MD) and no detection (ND), and the 
’’null” polarization and its error in gauss. 


# 

mid-HJD 

B, 

crB; 

Detect. 

N 

crN 

1 

2454391.559 

-5.7 

7.7 

MD 

-3.9 

7.7 

2 

2454392.719 

-26.9 

16.6 

ND 

37.2 

16.6 

3 

2454393.570 

-9.3 

5.5 

ND 

3.6 

5.5 

4 

2454394.491 

-0.3 

5.4 

ND 

6.3 

5.4 

5 

2454395.518 

18.1 

5.2 

MD 

12.4 

5.4 

6 

2454397.496 

25.5 

5.1 

MD 

3.9 

5.1 

7 

2454398.526 

4.7 

5.1 

MD 

-3.7 

5.1 

8 

2454762.644 

-15.1 

5.5 

MD 

-6.9 

5.5 

9 

2454763.645 

12.8 

6.6 

ND 

-5.4 

6.6 

10 

2454764.654 

28.0 

5.3 

MD 

3.6 

5.3 

11 

2454765.639 

32.8 

6.3 

DD 

11.1 

6.3 

12 

2454766.635 

10.9 

6.4 

MD 

12.7 

6.4 

13 

2455839.688 

-3.9 

6.5 

ND 

2.9 

6.5 

14 

2455840.670 

3.9 

6.5 

ND 

-6.7 

6.6 

15 

2455845.608 

21.3 

9.4 

ND 

4.8 

9.4 

16 

2455846.632 

-9.7 

11.0 

ND 

-1.8 

11.0 

17 

2455865.712 

12.4 

6.5 

ND 

-0.9 

6.6 

18 

2455873.557 

25.2 

13.5 

MD 

-3.6 

13.5 

19 

2455877.626 

13.6 

7.5 

ND 

-13.5 

7.5 

20 

2455878.565 

7.7 

7.5 

ND 

6.8 

7.5 

21 

2455890.673 

6.1 

8.7 

ND 

11.0 

8.7 

22 

2455891.660 

24.0 

7.4 

MD 

-6.7 

7.4 

23 

2455892.502 

6.0 

8.1 

ND 

-11.4 

8.1 

24 

2455895.667 

1.3 

6.9 

ND 

-8.0 

6.9 

25 

2455896.600 

-22.6 

22.3 

ND 

22.3 

22.5 

26 

2455910.477 

4.1 

13.8 

ND 

13.6 

13.8 

27 

2455935.555 

-7.2 

7.0 

MD 

4.2 

7.0 

28 

2455940.536 

5.9 

7.2 

ND 

-10.1 

7.2 

29 

2455941.539 

-3.1 

6.8 

MD 

4.5 

6.8 

30 

2455942.475 

-5.3 

8.0 

ND 

-4.4 

8.0 

31 

2455943.367 

4.8 

8.4 

ND 

4.3 

8.4 

32 

2455952.529 

25.2 

7.3 

MD 

-12.5 

7.3 

33 

2455953.431 

72.3 

59.1 

DD 

-60.4 

59.1 

34 

2455966.472 

51.0 

9.5 

MD 

-0.3 

9.5 

35 

2455967.402 

1.7 

8.5 

ND 

11.3 

8.4 

36 

2455968.343 

10.7 

16.7 

ND 

-13.7 

16.6 


23oct07 24oct07 






Fig. 6. Examples of LSD Stokes I (bottom), Stokes V (top), and 
null N (middle) profiles, normalized to Ic, for the secondary 
component ^ Ori Ah for a few nights of observations, computed 
from Narval data using only the 82 lines belonging to the sec¬ 
ondary component. The red line is a smoothed profile. 


culated the rate of detections in the 1000 models for each of the 
profiles depending on the field strength (see Fig. |^. 

We required a 90% detection rate to consider that the field 
would statistically be detected. This translates into an upper limit 
for the possible undetected dipolar field strength for each spec¬ 
trum, which varies between ~900 and ~2350 G (see Fig.|^. 

Since 36 spectra are at our disposal, statistics can be com¬ 
bined to extract a stricter upper limit taking into account that the 
field has not been detected in any of the 36 observations (see 
Neiner et al.|2015 1. The final upper limit derived from this com¬ 
bined probability for ^ Ori Ah for a 90% detection probability is 
~300 G (see thick line in Fig.0. 

7. Magnetic field configuration 

7.1. Rotational modulation 

We searched for a period of variation in the 36 longitudinal mag¬ 
netic field measurements of ^ Ori Aa with the clean-NG algo¬ 
rithm (see Gutierrez-Soto et al.|2009 1 . We obtained a frequency 
/ = 0.146421 c d ‘, which corresponds to a period of 6.829621 


days. This value is consistent with the period of ~1 days sug¬ 
gested by Bouret et al. ( 2008| l. Assuming that the magnetic field 
is a dipole with its axis inclined to the rotation axis, as is found 
in the vast majority of massive stars, this period corresponds to 
the rotation period of the star. 

We used this period and plotted the longitudinal magnetic 
field as a function of phase. For the data taken in 2007 and 2008, 
the phase-folded field measurements show a clear sinusoidal be¬ 
havior, as expected from a dipolar field model (see top panel 
of Fig. 1^. A dipolar fit to the data, that is, a sine fit of the 
form B{x) - Bq + Ba x sm(2ji{x + (pd)), resulted in Bq-6.9 G 
and Ba-19.2 G. A quadrupolar fit to the phase-folded data only 
shows an insignificant departure from the dipolar fit. 

However, the period of ~6.829 days does not to match the 
measurements collected in 2011 and 2012 very well (see middle 
panel of Fig. |^. None of the dipolar or quadrupolar fits to these 
data provide a reasonable match. A further search for a differ¬ 
ent period in these 2011-2012 data alone provided no significant 
result. 

The magnetic fields of main-sequence massive stars are of 
fossil origin. These fields are known to be stable over decades 
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Table 4. Longitudinal magnetic field measurements for the sec¬ 
ondary ^OriAb. The columns list the heliocentric Julian dates 
(HJD) for the middle of observation, the longitudinal magnetic 
field and its error in gauss, the detection status: no detection 
(ND) in all cases, and the ’’null” polarization and its error in 
gauss. 


# 

mid-HJD 

B, 

crB; 

Detect. N 

crN 

1 

2454391.559 

-62.9 

74.5 

ND 

46.5 

74.4 

2 

2454392.719 

-80.4 

132.7 

ND 

91.9 

132.8 

3 

2454393.570 

26.0 

54.8 

ND 

-60.6 

54.9 

4 

2454394.491 

31.0 

56.7 

ND 

24.1 

56.6 

5 

2454395.518 

86.7 

47.6 

ND 

101.4 

47.7 

6 

2454397.496 

-49.4 

45.5 

ND 

14.3 

45.9 

7 

2454398.526 

14.2 

46.0 

ND 

44.6 

46.1 

8 

2454762.644 

-14.2 

50.8 

ND 

73.0 

50.6 

9 

2454763.645 

-24.1 

65.4 

ND 

43.3 

65.5 

10 

2454764.654 

-40.1 

49.2 

ND 

80.0 

49.4 

11 

2454765.639 

-32.9 

73.5 

ND 

47.2 

73.6 

12 

2454766.635 

7.9 

60.5 

ND 

-26.5 

60.6 

13 

2455839.688 

-62.3 

53.2 

ND 

-30.5 

53.6 

14 

2455840.670 

-8.3 

58.3 

ND 

3.4 

58.2 

15 

2455845.608 

-19.1 

92.2 

ND 

-21.9 

91.1 

16 

2455846.632 

-126.4 

91.9 

ND 

-102.191.2 

17 

2455865.712 

35.7 

63.2 

ND 

173.4 

63.9 

18 

2455873.557 

141.5 

112.7 

ND 

-224.0113.0 

19 

2455877.626 

55.7 

60.0 

ND 

2.5 

60.0 

20 

2455878.565 

34.8 

68.3 

ND 

75.0 

68.7 

21 

2455890.673 

-105.9 

106.6 

ND 

-154.9107.5 

22 

2455891.660 

8.7 

82.6 

ND 

86.6 

82.8 

23 

2455892.502 

77.1 

84.3 

ND 

-63.7 

84.8 

24 

2455895.667 

-27.0 

75.2 

ND 

-63.8 

75.5 

25 

2455896.600 

-33.7 

289.4 

ND 

-113.7 296.7 

26 

2455910.477 

171.0 

183.7 

ND 

129.3 

183.4 

27 

2455935.555 

-82.4 

70.1 

ND 

-58.6 

69.9 

28 

2455940.536 

-12.3 

63.5 

ND 

-55.8 

63.7 

29 

2455941.539 

20.3 

67.3 

ND 

116.9 

67.8 

30 

2455942.475 

31.8 

70.6 

ND 

-97.2 

71.1 

31 

2455943.367 

-6.2 

78.8 

ND 

-27.7 

79.2 

32 

2455952.529 

-72.4 

71.0 

ND 

1.1 

71.3 

33 

2455953.431 

173.6 

478.5 

ND 

-152.6480.7 

34 

2455966.472 

52.9 

79.9 

ND 

-32.1 

79.8 

35 

2455967.402 

-57.8 

74.6 

ND 

62.8 

74.3 

36 

2455968.343 

-152.8 

186.5 

ND 

-92.9 

187.2 


and are only modulated by the rotation of the star. A change of 
period in the field modulation between the 2007-8 and 2011-12 
epochs is thus not expected in ^ Ori Aa. 

^ Ori Aa has a companion, therefore we investigated the pos¬ 
sibility that the magnetic field has been affected by the com¬ 
panion. Indeed, in 2011 and 2012, ^ Ori Ah was close to perias- 
tron, which means that the distance between the two stars was 
smaller than in 2007 and 2008. We calculated this distance to 
check whether some binary interactions might have occurred. 

To calculate the distance between the two components, we 


used the photometric distance of ^OriA, d=387 pc (Hummel 
|et al.||20l3| l. From [Hummel et ^ ( |201 S] !, we know the orbital 
parameters of the binary. The shortest distance between the two 


stars is - a- - b^- 23.8 mas, where a is the semi-major 
axis and b the semi-minor axis. From the distance of ( Ori A, we 
can compute rmin = sm(ff)/d in pc, where 0 is the parallax in 
radian. We obtained a distance of 81 R*, where R* is the radius 
of ( Ori Aa. 

The distance between the two stars at periastron therefore 
appears too large for interactions between the two stars to oc¬ 
cur. In addition, the binary system is still significantly eccentric 



Bpol (G) 

Fig. 7. Detection probability of a magnetic field in each spec¬ 
trum of the secondary component of ( Ori Ah (thin color lines) 
as a function of the magnetic polar field strength. The horizontal 
dashed line indicates the 90% detection probability, and the thick 
black curve (top left corner) shows the combined probability. 


(0.338, Hummel et al.||20T3] l even though f Ori Aa has already 
evolved into a supergiant. Tidal interactions have apparently not 
been able to circularize the system yet, which would confirm that 
these interactions are weak ( Zahn|2008l l. 

However, in addition to ^ Ori Aa and Ah, a third star f Ori B 
may also interfere with the ^ Ori A system. Correia et al. (2012 1 
showed that when a third component comes into play, tidal ef¬ 
fects combined with gravitational interactions may increase the 
eccentricity of ^OriA, which would otherwise have circular¬ 
ized. We thus cannot exclude that tidal interactions are stronger 
than they seem in the ( Ori A system. 


7.2. Field strength and geometrical configuration 


Assuming that the period detected in Sect. |7.1| is the rotation pe¬ 
riod of the star, we can determine the inclination angle of the star 
i by measuring v sin /. In massive stars, line broadening does not 
come from rotational broadening alone, but also from turbulence 
and stellar wind. This is particularly true for supergiant stars. 


Based on the synthetic spectra calculated in Sect. |4.1 [ and the 
lines identified to belong to only one of the two components, 
we determined the broadening needed in the synthetic spectra 
to fit the observations. For f Ori Aa, a broadening of 230 km s ' 
was necessary to provide a good fit to the observed lines, while 
for ^ Ori Ah we needed 100 km s '. These broadening values are 
upper limits of the v sin i values because they include all physical 
processes that broaden the lines. In fact, with a period of 6.829 
days and a radius of 20 Rq as given by Hummel et al. P013| l, the 
maximum possible v sin i for ( Ori Aa is 148 kms L 


In addition, [Bouret et al.| ( |2008l l determined vsini through a 
Fourier transform of the average of the 5801 and 5812 A CIV 
and 5592 A OIII line profiles. They found a vsin/ of 110+10 
km s'. From our disentangling of the spectra, we know that the 
two CIV lines originate from ( Ori Aa, but the OIII 5590 A line 
is partly (~ 10%) polluted by ^ Ori Ah. As a consequence, we ap¬ 
plied the Fourier transform method to the LSD Stokes I profiles 
we calculated from the lines that only originate from ^ Ori Aa. 
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Fig. 8. Rotational modulation of the longitudinal magnetic field 
of ^OriAa for the observations taken in 2007-2008 (top) and 
2011-2012 (center). The black line corresponds to the best 
dipolar fit, while the dashed blue line corresponds to the best 
quadrupolar ht. The bottom panel compares the fit of the dipole 
and the quadrupole obtained from the 2007-2008 data with the 
observations obtained in 2011-2012 (see Sect. |7.3[ l. The data for 
which the Stokes V model matches the observed LSD V profiles 
are shown in black, while the data for which the Stokes V model 
does not match are in red. 


We obtained v sin /= 140 km s'. However, it is known that v sin i 
values determined from LSD prohles might be overestimated. 

Finally, taking macroturbulence into account but not binarity, 
foe example, Simon-Diaz & Herrero (|2014| l found that v sin i for 
( Ori A is between 102 and 127 km s ', depending on the method 
they used. 

In the following, we thus consider that v sin i is between 100 
kms 'and 148 kms 'for ^Ori Aa. In addition, we adopt the ra¬ 


dius of 20 Rq given by [Hummel et al.| ( |2013j l and the rotation 
period of 6.829 d. Using vsini= [100-148] km s ', we obtain 
i ~ [42-87]°. 

Using the dipolar fit to the 2007-2008 longitudinal field 
measurements and the inclination angle i, we can deduce the 
obliquity angle (3 of the magnetic field with respect to the ro¬ 
tation axis. To this aim, we used the formula r - Bmin/Smax = 
cos(/ - yS)/ cos(! + P) ( |Shore|1987 1 . The dipolar fit of the longi¬ 
tudinal field values gives r=0.47. With i ~ [42 - 87]°, we obtain 
P~ [71-8]°. 

In addition, from the dipolar ht to the longitudinal held 
values and the angles i and p determined above, we can es¬ 
timate the polar held strength with the formula Bq ± Bg - 
0.296 X Bpoi cos(jS + i), where the limb-darkening coefhcient is 
assumed to be 0.4 (see Borra & Landstreet 1980 1 . We found Bpoi 
= [110 + 5-524 + 65] (j. Ihe dipolar magnetic held that we hnd 


is thus higher than the one found by Bouret et al. (2008 1 . 

In 2011-2012, the maximum measured B; is 51 G and the 
minimum polar held strength is thus Bpoi > 33Bi^max = 168 + 
33 G. This value is compatible with the range derived from the 
2007-2008 data. 


7.3. Stokes V modeling 

Since the Bi data taken in 2007-2008 point toward the presence 
of a dipole held, we used an oblique rotator model to ht the LSD 
Stokes V and I prohles. 

We used Gaussian local intensity prohles with a width calcu¬ 
lated according to the resolving power of N arval and a macrotur - 
bulence value of 100 km s ' determined by Bouret et al. (2008 i. 
We ht the observed LSD I prohles by Gaussian prohles to de¬ 
termine the depth, v sin i and radial velocity of the intensity pro- 
hle. We used the weighted mean Lande factor and wavelength 
derived from the LSD mask applied to the Narval observations 
and the rotation period of 6.829 days. The ht includes hve pa¬ 
rameters: i, p, Bpoi, a phase shift (p, and a possible off-centering 
distance dj of the dipole with respect to the center of the star 
(dj=0 for a centered dipole and d^=l if the center of the dipole 
is at the surface of the star). 


8. Magnetospheres 

8.1. Magnetospheric parameters 

We calculated a grid of V prohles for each phase of observation 
by varying the hve parameters mentioned above and applied a 
minimization to obtain the best ht of all observations simulta¬ 
neously. More details of the modeling technique can be found in 
Alecian et al.|( 20081. The parameters of the best ht are i=79.89°. 


P-21.5°, 0=0.68, Bpoi=142.2 G and d^=0.0. The values for the 
angles i and p are within the error boxes derived in Sect. 7.2 and 
the value for the polar held strength Bpoi hts the Bi results well. 
Moreover, the best ht is obtained for drf=0, which conhrms that 
no quadrupolar component is found. 

The 36 Stokes V prohles for this best ht are shown in Fig. 
overplotted on the observations. As expected, for the nights in 
2007-2008, the model hts the observations well. For some nights 
in 2011-2012, the observations are too noisy to see whether the 
model hts well. Considering the nights in 2011-2012 for which 
the S/N is sufficient, the model hts some the observations but 
not all. For those nights for which the model htted well the ob¬ 
servations, we compared the values of the longitudinal magnetic 
held B; to the dipolar ht obtained for the B/ measurements of 
2007-2008 (see bottom panel of Fig. [^. The 2011-2012 data that 
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Fig. 9. Best dipolar model fit (green) of the observed Stokes V profiles (black). The green numbers correspond to the rotational 
phase. The very noisy observation at phase 0.793 is not shown. 


match the Stokes V models also match the Bi dipolar fit curve. 
Therefore, it seems that at least part of the 2011-2012 data show 
the same rotational modulation and dipole field as in 2007-2008. 
Only part of the 2011-2012 dataset does not match the rest of the 
observations. 


With the polar magnetic field strength Bpd = 142.2 G de¬ 
termined with to the Stokes V model, we calculated the wind 
confinement parameter /;» which characterizes the ability of the 
magnetic field to confine the wind particles into a magneto¬ 
sphere (ud-Doula & Owocki 2002 1 . If < 1, ^OriAa is lo¬ 
cated in the weakly magnetized winds region of the magnetic 


confinement-rotation diagram and it does not have a magneto¬ 
sphere. However, for ?/» > 1, wind material is channeled along 
magnetic field lines toward the magnetic equator and ^ Ori Aa 
hosts a magnetosphere. 


To calculate rj,, we first used the fiducial mass-loss rate 
= 1.4 X IO^^Mq and the termin al speed Voo = 
2100 km s“* determined by Bouret et al. (2008 1 . They measured 


the mass-loss rate from the emission of Ha and used archival 
International Ultraviolet Explorer (lUE) spectra to measure the 
wind terminal velocity from the blueward extension of the strong 
UV P Cygni profile. We obtain 77 * = 0.9. 
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Fig. 10. Dynamic plot of each individual archival Hq, spectrum 
in phase with the rotation period. 


We then recalculated 77 * but this time using the mass-loss rate 
of M - 3.4 X IO^^Mq yr~^ and Voo = 1850 kms“^ determined 


by Cohen et al. (2014 1 . This gives 77 , = 4.2. 

A magnetosphere can only exist below the Alfven ra¬ 
dius Ra, which is proportional to t/*, with Ra - T]],^^Rf For 
^OriAa, using the two above determinations of rj^, Ra - 
[0.98 - 1.43]/?*. Moreover, the magnetosphere can be centrifu- 
gally supported above the corotation Kepler radius Rk- Rk - 
(27r/?t/Prot thus for ^Ori Aa Rk - 2.8/?*. Since 

Rk I Ra, no centrifugally supported magnetosphere can exist. 

Therefore, ( Ori Aa is either in the weakly magnetized winds 
region of the magnetic confinement-rotation diagram, meaning 
that ( Ori Aa does not have a magnetosphere ( 77 * < 1), or it hosts 
a dynamical magnetosphere (1 < 77 * < 4.2). 


8.2. Ha variations 

The H„ line shows significant variability in emission and ab¬ 
sorption. For stars that have a magnetosphere, we expect magne- 
tospheric emission at Hq,, which varies with the rotation period 
(see e.g. |Grunhut|2015| l. 

To check whether there is a signature of the presence of a 
magnetosphere around f Ori Aa, we studied the variation of its 
Hq line in the archival spectra (see Sect. \2.2\ . We confirm that 
the emission in H„ does indeed vary. While most of the vari¬ 
ations are problably related to variations in the stellar wind of 
the supergiant, the signature of a weak rotationally modulated 
dynamical magnetosphere is observed in Hq (see Fig.[T 0 |l. 

The ratio \ogRAlRK gives a measure of the volume of the 
magnetosphere. For ^Ori Aa, logRAlRK is very small (< 0.06), 
and it is thus not surprising that Hq only weakly reflects mag¬ 
netic confinement. 


9. Discussion and conciusions 


Based on archival spectrocopic data and Narval spectropolari- 
metric data, we confirm the presence of a magnetic held in 
the massive star (Ori A, as initally suggested by [Bouret et al. 
( 2008| l. However, Bouret et al. ( 2008| l ignored that ^OriA is a 
binary star, which was subsequently shown by [Hummel et al.| 
( 2013| l with interferometry. 

We disentangle the spectra and could thus show that the 
primary O supergiant component ( Ori Aa is the magnetic star, 
while the secondary ^ Ori Ah is not magnetic at the achieved de¬ 
tection level. (, Ori Aa is the only magnetic O supergiant known 
as of today. 

The magnetic held of ^OriAa is a typical oblique dipole 
held, similar to those observed in main-sequence massive stars. 
From Stokes modeling, the polar magnetic held strength Bpoi of 
^ Ori Aa is found to be about 140 G. If we assume held conser¬ 
vation during the evolution of ^ Ori Aa because the stellar radius 
increased from ~10 to ~20 Rq, the surface magnetic polar held 
strength decreased by a factor ~4. This implies that the polar 
held strength of ^ Ori Aa when it was on the main sequence was 
about 600 G. This is similar to what is observed in other main- 
sequence magnetic O stars. 

The current held strength and rotation rate of (Ori Aa are 
weak, with respect to the wind energy, for the star to be able 
to host a centrifugally supported magnetosphere. However, it 
seems to host a dynamical magnetosphere. All other ten known 
magnetic O stars host dynamical magnetospheres, except for the 
complicated system of Plaskett’s star, which has a very strong 
magnetic held and hosts a centrifugally supported magneto¬ 
sphere (see jGrunhut et al.|2013| l. However, these other magnetic 
O stars are not supergiants. 

Although ^ Ori A is one of the brightest O star in the X- 
ray domain, [Cohen et al.j ( [2014[ l found that it resembles a non¬ 
magnetic star, with no evidence for magnetic activity in the X- 
ray domain and a spherical wind. This probably results from the 
weakness of the magnetosphere around ^ Ori Aa. 

The rotation period of ^'OriAa, - 6.829 d, was deter¬ 
mined from the variations of the longitudinal magnetic held. 
This period is clearly seen in the data obtained in 2007-2008, 
but only part of the spectropolarimetric measurements obtained 
in 2011 and 2012 seem to follow that rotational modulation. The 
reason for the lack of periodicity for part of the magnetic mea¬ 
surements of 2011-2012 was not identified. Although passage at 
the binary periastron occured between 2008 and 2011 , the dis¬ 
tance between the two companions seems too large for the com¬ 
panion to have perturbed the magnetic held of the primary star, 
unless it is f Ori B which has maintained the two components of 
^ Ori A at a distance (see Sect. HD- 

^ Ori A therefore remains an interesting star that needs to be 
studied further. More spectropolarimetric observations should be 
collected at appropriate orbital phases to allow for a more accu¬ 
rate spectral disentangling, this would allow obtaining stronger 
constraints on the magnetic held strength and configuration, 
studying the held as a function of orbital phase, and understand¬ 
ing the magnetic held perturbations that seem to have occurred 
during the observations in 2011 - 2012 . 
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Fig. 11. LSD Stokes I profiles (bottom) computed from the disentangled spectroscopic data, Stokes V (top) and null N (middle) 
profiles, normalized to Ic, from the Narval data, for ^ Ori A. The red line is a smoothed profile. 























































































































